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Abstract. We present a measurement of the charge asymmetry in top quark pair production using an integrated 
luminosity of 1.09 fhr 1 collected with the CMS detector. Top quark pairs with a signature of one electron or 
muon and four or more jets, at least one of them b tagged, are selected. At the LHC a small charge asymmetry 
in the rapidity distributions of top and antitop quarks is predicted. Therein slightly broader rapidity distributions 
for top quarks are expected, while antitop quarks are produced more centrally and possess narrower rapidity 
distributions. We determine the charge asymmetry based on two different sensitive variables and the results are 
compared with the most precise standard model theory predictions using a dedicated unfolding technique. 



1 Introduction 

In the Standard Model (SM), a small charge asymmetry 
in tt production via quark-antiquark annihilation is pre- 
dicted. It can be explained via the interference between the 
Born and the box diagram, and between initial- and final- 
state radiation, which link the flight directions of the top 
quarks and antiquarks to the directions of the initial quarks 
and antiquarks, respectively 111 I2I . At the Tevatron proton- 
antiproton collider this effect leads to a forward-backward 
asymmetry which is visible in the difference in rapidity 
y of top quarks and antiquarks. Recent measurements by 
the CDF and D0 collaborations [fFfl measure asymme- 
tries which are about two standard deviations larger than 
the predicted SM value of about 0.08 0. The CDF collab- 
oration reports an even larger asymmetry compared to the 
SM prediction at high tt invariant masses |6 1. These results 
have led to speculations that the large asymmetry might be 
generated by additional axial couplings of the gluon or by 
heavy particles with unequal vector and axial-vector cou- 
plings to top quarks and antiquarks. 

In proton-proton collisions at the LHC the charge asym- 
metry leads to an central-peripheral asymmetry, which can 
be observed through the difference in the absolute values 
of the pseudorapidities of top quarks and antiquarks A\r\\ — 
M - We also measure the charge asymmetry in the ob- 
servable used by the Tevatron experiments multiplied by a 
factor that accounts for the boost of the tt system, yielding 
4y 2 = (y t - m) • iVt + m) = (yf - motivated in 0. Using 
either of the two variables, the charge asymmetry can be 
defined as 



where and N are the number of events with positive 
and negative values in the sensitive variable, respectively. 
The SM prediction for^%| isA£(theory) = 0.0136+0.0008 
and for Ay 2 it is A^(theory) = 0.0115 ± 0.0006 0, 
where the uncertainties arise from scale variations and par- 
ticle density function uncertainties. These values are much 
smaller than the SM predictions for the Tevatron, because 
of the large fraction of gluon-gluon induced tt production. 



2 The CMS Detector 

The central feature of the Compact Muon Solenoid (CMS) 
apparatus is a superconducting solenoid, of 6 m internal 
diameter, providing a field of 3.8 T. Within the field vol- 
ume are the silicon pixel and strip tracker, the crystal elec- 
tromagnetic calorimeter (ECAL) and the brass/scintillator 
hadron calorimeter (HCAL). Muons are measured in gas- 
ionization detectors embedded in the steel return yoke. In 
addition to the barrel and endcap detectors, CMS has ex- 
tensive forward calorimetry. A much more detailed descrip- 
tion of CMS can be found elsewhere (8). 



3 Selection of Events 

The measurement is performed in the electron+jets and 
muon+jets decay channels, using proton-proton collision 
events at a center of mass energy of 7 TeV recorded with 
the CMS detector up to July 201 1 . The amount of data cor- 
responds to an integrated luminosity of 1.09 ± 0.07 fb _1 . 
Only events selected by special trigger algorithms which 
search for electrons or muons together with at least three 
jets with at least 30GeV/c of transverse momentum are 
used. In addition, the presence of a good primary vertex 
(PV) is required. 

In the electron+jets channel electron candidates are re- 
quired to possess a transverse energy larger than 30 GeV 
and to lie in the region \tj\ < 2.5, whereas in the muon+jets 
channel muon candidates must have a transverse momen- 
tum greater than 20GeV/c and must lie within the muon 
trigger acceptance (\r]\ < 2.1). Both lepton candidates have 
to be isolated and must pass additional quality criteria. Jets 
are required to have a corrected jet energy p T > 30GeV/c 
and \t]\ < 2.4. 

Only events with exactly one isolated charged lepton 
and at least four jets, where one of them is tagged as a 
b jet, are selected. In both channels events with additional 
more loosely defined charged leptons are discarded. 

Applying the event selections to the recorded dataset 
of 1.09 ft)" 1 in total 12757 events are selected, 5665 in the 
electron+jets channel and 7092 in the muon+jets channel. 
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Fig. 1. Reconstructed A\rj\ distribution for the combined 
lepton+jets channel. 

4 Background Estimation 

The measurement of the charge asymmetry requires a pre- 
cise knowledge of the amount if background events in the 
analyzed dataset. Therefore we apply a data-driven back- 
ground estimation using two kinematic variables, the miss- 
ing transverse energy and M3, which is the invariant mass 
of the three jets with the largest vectorially summed trans- 
verse momentum. The selected datasets in the electron+jets 
and muon+jets channel are separated each in two subsets, 
events with < 40GeV and events with pj > 40GeV. 
We fit the distribution in the low-pj dataset and the 
M3 distribution in the high-^j dataset simultaneously. The 
background estimation is performed individually for the 
two lepton types. W+jets contributions from W + and W~ 
are estimated individually. In order to simulate the QCD 
background properly a QCD model from data is used by 
inverting the cut on the relative isolation. The fits are per- 
formed as binned likelihood fits using the theta frame- 
work [9 1 yielding a signal purity of about 80 %. 



5 Reconstruction of tt Pairs 

In order to fully reconstruct the top quarks in the com- 
bined lepton+jets decay channel, the selected physical ob- 
jects are consecutively assigned to final state particles. The 
assignment is not unambiguous and we define a criterion 
iff based only on reconstructible quantities in order to give 
each possible hypothesis a probability to be the best hy- 
pothesis for the event. For t// the three linearly decorre- 
lated masses nt\, tn%, and arising from the masses of 
the reconstructed top quarks and the hadronically decay- 
ing W boson, as well as the b tagger discriminator values 
for the jets assigned to the two b quarks and to the two 
light quarks are taken into account. The hypothesis with 
the biggest value of if/ is then chosen in each event for fur- 
ther consideration. 



6 Measurement of the tt Charge 
Asymmetry 

In figure[T]the distribution of A\r]\ in the lepton+jets chan- 
nel obtained from the reconstructed top and antitop quarks 
is shown. Both reconstructed distributions (A\rj\ and Ay 2 ) 
can be used to calculate an uncorrected charge asymmetry 



A^. AW by simply using the definition in equation|l|yielding 
A c aw (tj) = -0.004 ± 0.009 and A R c AW (y) = -0.004 ± 0.009. 

These values are not directly comparable with any the- 
oretically motivated prediction, since several effects bias 
the measurement at this stage. First of all, 20% of events 
used to measure A C AW arise from background processes. 
Therefore we subtract the predicted amount of background 
events from the measured distributions. In this step a Gaus- 
sian error propagation is performed taking the measured 
uncertainties on the background rates as input. 

An even larger effect is due to the event selection ef- 
ficiency and to the imperfect reconstruction method. De- 
pending on the true value of A\rj\ (or Ay 2 ) the probability 
for an event to survive the event selection varies, which 
leads to distortions in the reconstructed distributions. Fur- 
thermore, in many events not exactly the best possible jet 
parton assignment and neutrino reconstruction is chosen. 
But even when the best possible reconstruction hypothe- 
sis is chosen, this does not guarantee that the top quark 
four-momenta are reconstructed correctly. Due to the fi- 
nite energy resolution of the calorimeters and the jet re- 
construction, and due to the possibility that jets from the 
tt decay can lie outside the detector acceptance, a perfect 
reconstruction is a priori not possible. In order to investi- 
gate all these effects one can analyze simulated tt events, 
where the four momenta of the generated top quarks are 
compared to the four vectors reconstructed from measured 
objects like jets and leptons. From these studies one finds 
that the distortion of the distributions can be factorized into 
effects due to the event selection efficiency and due to the 
event reconstruction procedure. 

To correct for these effects, we apply a regularized un- 
folding procedure IfTUll . The unfolding algorithm corrects 
the measured spectrum for migration and efficiency effects 
by applying a generalized matrix inversion method. The 
measured spectrum is divided into 12 bins, where the dif- 
ferent widths of the bins have been chosen such that they 
contain approximately equal numbers of events. For the 
unfolded spectrum six bins have been used. To regularize 
the problem and to avoid unphysical fluctuations two ad- 
ditional terms, a regularization term and a normalization 
term, are used 01 1 II 21 . 

The performance of the unfolding algorithm is tested in 
50,000 pseudo experiments. For each pseudo experiment a 
pseudo dataset is created and unfolded with the method 
described above. In each pseudo experiment the unfolded 
spectrum is compared with the generated distribution and 
a very good agreement is found putting the analysis on 
solid ground. Also with pseudo experiments in which we 
re-weight the events in order to generate a non-zero asym- 
metry between -0.21 and +0.23 a linear behavior of the 
average fit results as a function of the true values is found. 
While for the A\rj\ variable the agreement is almost perfect, 
for Ay 2 a small bias in the order of 0. 1 % has been observed 
for which the measured asymmetry has to be corrected by 
a factor of 0.94. 

The measurement of the charge asymmetry Ac might 
be affected by several systematic sources uncertainties. Only 
systematic uncertainties influencing the direction of the re- 
constructed top quark momenta can change the value of 
the reconstructed charge asymmetry. The overall selection 
efficiency and acceptance will not change the measured 
values. For each source of systematic pseudo datasets are 
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Fig. 2. Unfolded A\t]\ normalized spectrum. The NLO prediction 
is based on the calculation of Ref. (3). 
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Fig. 3. Background-subtracted asymmetries for A\rj\ as function 
of the reconstructed tt invariant masses. 

drawn from systematically shifted samples, to which the 
unfolding method with the standard templates are applied. 
The largest systematic uncertainties arise from the varia- 
tion of the Q 2 scale and matching threshold and from the 
variation of initial- and final-state-radiation in the used tt 
signal simulated sample. 

We apply the described unfolding procedure to the mea- 
sured A\rj\ distribution as well as to the distribution of the 
second variable, Ay 2 . Figure [2] shows the unfolded spec- 
trum for z/|?7| used for computing the asymmetries together 
with the SM prediction at NLO. In the unfolded A\j]\ distri- 
bution we measure an asymmetry of 

Al = -0.016 + 0.030 (stat.)!™|°(syst.) , (2) 

while in Ay 2 we measure an unfolded and corrected (di- 
vided by 0.94) asymmetry of 

A y c = -0.013 + 0.026 (stat.)_^(syst.) . (3) 

Both measured values are within the uncertainties in 
agreement with the theory predictions. 

One can also measure the background subtracted asym- 
metry as a function of the reconstructed invariant mass of 
the tt system to investigate whether one can see a depen- 
dence of the asymmetry on ma. Figure [3] shows the result 
for the A\rj\ variable, where no increase of the asymmetry 
for increasing m t j can be seen. The same behavior is found 
for Ay 2 . However, these studies allow only for a qualita- 
tive statement, while for a quantitative statement a proper 
simultaneous unfolding in the sensitive variable as well as 
in the invariant mass of the tt system has to be performed. 



7 Conclusion 

The analysis iTPJI presented here provides the measurement 
of the charge asymmetry in tt events using a dataset corre- 
sponding to an integrated luminosity of 1.09 fb -1 . A se- 
lection of tt events has been performed in the lepton+jets 
decay channel and the amount of background events in the 
selected dataset has been estimated via a binned likelihood 
fit. Thereafter the four-momenta of the top quarks have 
been reconstructed and corrected using a regularized un- 
folding method. As opposed to the Tevatron results, we 
measure a small negative asymmetry which is within the 
uncertainties still very well compatible with the SM pre- 
dictions and shows no tendency to large deviations from 
the prediction. Also the background subtracted asymmetry 
as a function of the invariant mass of the tt system shows 
no tendency to larger values for large invariant masses. For 
a quantitative statement of this behavior a proper simul- 
taneous unfolding in both, the sensitive variable and the 
invariant mass of the tt system, will be performed in the 
future in order to confirm or to refute the Tevatron find- 
ings. 
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